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An all optical soliton diode (AOSD) is proposed based on a sandwich nonlinear Bragg-grating 
structure: a linearly chirped Bragg-grating linked to a uniform Bragg-grating and again to a chirped 
Bragg-grating. The nonreciprocity is achieved by introducing two spatially asymmetric chirped 
Bragg-gratings with optical nonlinearity. High transmission ratio up to 150 is obtained when launch¬ 
ing a picosecond Gaussian pulse into the setting. We find that such pulses in a form of solitons 
propagate at a rather small velocity (< 0.03c) when pulse wavelength is selected in the vicinity of 
photonic bandgap. 

OCIS code: (060.3735) Fiber Bragg gratings;(190.5530) Pulse propagation and temporal solitons; 
(230.4320) Nonlinear optical devices. 


Optical diode, a nonreciprocal device that transmits 
light only in one direction and blocks that in the opposite 
direction, is a key enabler for optical information process¬ 
ing. To realize the optical diode, breaking the symmetry 
of wave propagation is necessary but it is still a challenge 
problem. Many schemes have been proposed to solve 
such a problem. Conventional methods were based on the 
Faraday effect in magneto-optical crystals [l|, but such 
magneto-optical component hinders the applications in 
integrated photonic systems. Since then, the research on 
optical diodes became very intense, and many compact 
schemes were proposed, using opto-acoustic effect (2j, sec¬ 
ond harmonic generation [3|, absorbing multilayer system 
Q, asymmetric nonlinear absorption [5], ring-resonator 
structures 6, 7], Fano effect [8, 9], micro-optomechanical 
devices with flexible Fabry-perot [10[ , left-hand materials 
mi and photonic crystal structures nng. While ma¬ 
jorities of these work concentrated on high nonreciprocal 
transmission ratio [7, 15], the shape of the optical field 
after propagating through the devices was neglected. 

In addition to those attempts, an important candidate 
for the generation of optical diodes is based on periodic 
grating structures due to their tunability and cascadabil- 
ity. Using the edge of photonic bandgap mi , as well 
as the consideration of optical nonlinearity, the optical 
diode was first put forward in a one-dimensional pho¬ 
tonic material with spatial gradation in the refractive 
index. Later, an all-optical diode was demonstrated in 
a periodically LiNbCU channel waveguide using quasi¬ 
phase matching [18]. Non-reciprocal transmission of cir¬ 
cular polarized light at the photonic-bandgap regions 
based on electro-tunable liquid-crystal hetero-junctions 
was demonstrated experimentally and theoretically [19] . 

Similarly to 0, we propose a sandwich structure 
based on nonlinear Bragg gratings: a linearly chirped 
Bragg-grating (CBG) [20j linked through a uniform 
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FIG. 1: (Color online) An AOSD sketch constructed by a 
sandwich structure: CBG+UBG+CBG. 


Bragg-graitng (UBG), and again to a CBG. Using this 
setting, a hight-Q and slow-light all-optical soliton diode 
(AOSD) is achieved. It was demonstrated thaL by intro¬ 
ducing a CBG cascaded a UBG structure [2l], [22[, strong 
retardation of optical pulses, as well as the creation of 
standing stable pulses were realized. These pulses can 
preserve its shape while propagating along the gratings, 
as a results of the balance between optical induced non¬ 
linearity and pulse dispersion, which can be termed as 
a BG soliton pdl l3lj. In our setting, it is possible to 
break the reciprocity by introducing such two spatially 
asymmetric CBG structures, linked via a UBG segment. 
The paper is organized as following, we firstly present 
the standard nonlinear coupled-mode theory to discuss 
the possibility for the realization of the nonreciprocity of 
the device. Then we apply such nonreciprocity to realize 
an AOSD. Performance of such diodes is systematically 
discussed. Conclusions are made in the end of this paper. 

Fig. m shows the sketch of the sandwich struc¬ 
ture. The total length of the setting is L, and the 
lengths of two CBG segments are rq and rq, respec¬ 
tively. The linear refractive index (2l|,[22|: n[z ) = no[l + 
2 Ancos( 27 rz(l+C±z)/ Ao)] for the left CBG segment (i.e., 
z G [0,ri]); n(z) = no[l + 2An cos(2^(1 + CAi/Ao)] 
for the UBG segment (i.e., z G (ri,L — 7*2)); n(z) = 
no[l + 2Ancos(27rz[l + C 2 (I/ — z)]/A 0 )] for the right CBG 
segment (i.e., z G [L — r 2 ,L]). Here, z is the propagation 
distance. C\ and C2 are the chirp coefficients of these 
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FIG. 2: (Color online) Typical examples of pulse propagation 
for two different length of CBG, n =0.1 cm (a) and 77 — 0.6 
cm (c), with the same PI, I = 4.24 GW/cm 2 . (b,d)Pf and 
Pb versus z, corresponding to (a,c) respectively. In (b,d), the 
gray areas are the CBG regions and the blue solid curves are 
the Pf(z)/Ptotaij while the red dash curves are Pb(z)/Ptotai- 
Here Ptotai = /^(IP/I 2 + \Eb\ 2 )dt. and L is fixed to 2 cm. 


two CBG segments, respectively. Ao is the BG period at 
the input and output edges, no is the average refractive 
index, and An is the amplitude of the refractive-index 
modulation. 

The BG solitons propagating through the setting are 
well described by the following nonlinear coupled-mode 
equations (CMEs) [H 12, |H| 


+ + 5{z)E f' b + K{z)EbJ 

+7(l- E ’/,b| 2 + 2\E bt f\ 2 )E f) b = 0, 


( 1 ) 


where t is time, Ef and E & are amplitudes of the forward¬ 
traveling and backward-traveling waves of the soliton. 
v g = c/tiq is the linear group velocity in the material 
of which the BG is fabricated, with c being the light 
speed in the vacuum, and 7 = ri 2 Uj/c is the nonlinearity 
coefficient, with uj being the carrier frequency and ri 2 the 
Kerr coefficient. The coupling strength between forward 
and backward waves is k = 7rAn/Ao, while S(z), which 
represents the local wavenumber detuning, is 


5(z) 


(6 0 - 2ttCiz/Aq 

l So ~ AS 

(7 - 27 tC 2 (L - z)/ A, 


z € [0, »*i), 

z €[r 1 ,L-r 2 ], ( 2 ) 

z e (L- r 2 , L\, 


where So = 27mo/A — tt/Ao is the detuning at the left and 
right edges of the CBG segment, and AS = 27rCi77/A 0 = 
27rC2r2/Ao, which is a constant, is the detuning differ¬ 
ence between the UBG segment and the left/right edge 
of the setting. 

Simulations of Eq. © were performed by means of the 
fourth-order Runge-Kutta method, with physical param¬ 
eters of the silicon, whose self-focusing Kerr coefficient is 
4.5 x 10 “ 14 cm 2 /W [33j], and the average refractive index 
is no = 3.42. The refractive index modulation depth is 


taken as An = 0.006. To avoid the two-photon absorp¬ 
tion [33j], the carrier wavelength of 1798.4 nm is selected, 
and the BG period at the input edge of the sample is 
fixed to Ao = 263.16 nm, giving rise to £0 = 106.2 cm -1 . 
At the input edges of the sample, a picosecond Gaussian 
pulse, with a temporal width 3.1 ps and a tunable peak 
intensity (PI), namely 7, is coupled into the system. 

As mentioned [22], the introduction of CBG segment 
decreases the pulse propagating group velocity, thus af¬ 
fecting the optical transmission. To see this clearly, we 
fix f 2 =0 and launch the pulse from the left edge of 
the sample (i.e. side A). Moreover, we fix AS = 7.31 
cm -1 , so that the considered wavelength is located at 
the edge of the photonic bandgap (here, the chirp coef¬ 
ficient, C\ = A5Ao/27rri, is a function of 7 * 1 ). Typical 
examples of the pulse transmission and reflection with 
two different values of r\ are displayed in Fig. [2] (a,c), 
under identical incident intensity. From Fig. ©a), it 
is observed that the pulses maintain its shape and pass 
through the setting. It indicates that strong optical non¬ 
linearity is induced, overcoming the pulse dispersion near 
the edge of bandgap. As a result, a BG soliton is gen¬ 
erated propagating in UBG segment. Moreover, as ex¬ 
pected, the pulse cannot penetrate the setting with a 
larger CBG length, see Fig. ©<0, eventually bounced 
back to the incident edge. These assertions can be ver¬ 
ified by the comparison between these two characters: 
Pb,i( z ) = I-oo l-^b.fO, t)\ 2 dt, shown in Fig. HKb,d), corre- 
sponding to Fig. ©a,c). From these panels, it is seen that 
Pf transfers to Pb during propagating inside the CBG 
segment. In panel. ©b), corresponding to a smaller 7 * 1 , 
Pf — Pb > 0, and becomes constant in UBG segment. Ac¬ 
cording to Up U i se = v g ( 1 - /)/(1 + /) ( where / = P b /Pf 
J22[), we have U pu ise > 0. That means the pulse keeps 
propagating in the forward direction through the setting. 
While in panel. ©d), corresponding to a larger 7 * 1 , it is 
observed that majority of injected power is converted into 
backward-wave before the pulse reaches the interface. As 
a result, the pulse bounces back. 

To identify the relationship between CBG length and 
the optical transmission, the transmissivity is investi¬ 
gated for two selected r \, with the following expression 


_ f[\E f (z = L,t)\ 2 + \E b (z = L,t)\ 2 ]dt 
71 J[\E f (z = 0,t)\ 2 + \E b (z = 0,t)\ 2 }dt 


x 100%. (3) 


As an example, Fig. [3fa) displays 7 , as a function of in¬ 
cident PI for 7*1 =0.1 (red solid curve) and 0.6 cm (black 
dash curve), respectively. From this figure, a threshold 
value 7 t h is found for a fixed 77 (when I < 7 t h, 7 ~ 0 ; 
while I > 7 t h, 7 rises significatly). In this case, we find: 
7th =4 GW/cm 2 for 77 = 0 . 1 cm; 7 t h =4.38 GW/cm 2 for 
77 = 0.6 cm. Note that when 7 is selected near 7 t h, a 
nearly stopped soliton propagating in UBG segment can 
be achieved, see Ref. j2l|. Moreover, the relation be¬ 
tween 7 t h and 77 is also investigated, see Fig. [3fb). As 
expected, because of the strong pulse retardation in CBG 
region, the increase of 77 also leads to an increase of the 
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FIG. 3: (Color online) (a) 77 versus /, the PI of the input 
pulse, for different n: 0.1 cm (red solid curve) and 0.6 cm 
(black dash curve), (b) 7 t h vs. r 1 . (c) Factor Q of the AOSD, 
as a function of incident intensities, (d) Velocity of the BG 
soliton, namely, V s , versus PI, which is roughly estimated by 
L/r, where r is the time when the soliton propagates from 
the left end to the right end (i.e. A B). Here, in all the 
panels, we fix AS = 7.31 cm -1 and L — 2 cm. 



FIG. 4: (Color online)A typical example of an isolator for the 
setting in Fig. |T] The PI of the pulse is I — 4.35 GW/cm 2 . 
(a) The pulse injected from side A (n = 0.1 cm), which can 
transmit the setting. For this case, V s ~ 2 cm/2.43 ns= 
0.0274c. (b) The pulse injected from side B (7*2 = 0.6 cm), 
which is bounced back by the setting. Here, we fix AS = 7.31 
cm -1 and L — 2 cm. 


ith. 

It is interesting to note that, see Fig. [3(a), the non¬ 
reciprocity can be achieved as long as another CBG 
segment is added to the right side, namely the swand- 
wich structure in Fig. 1. For example, given rq = 0.1 
cm, 7*2 = 0.6 cm, an all-optical diode for an unidirec¬ 
tional transmission of soliton could be realized: within 
4.0 GW/cm 2 < / <4.38 GW/cm 2 , when the pulse is in¬ 
jected into the setting from the left side (namely, side 
A), it can evolute to a soliton and penetrate the setting 
with high transmission; while the pulse incident from the 
right-side end (namely, side B), it is nearly completely 
blocked. 

To identify the quality of the AOSD, we introduce a Q 
factor of this AOSD as Q = tib/vauM, where t)a is the 
transmissivity of the soliton from B —>> A, while t]b is the 
transmissivity from A —>> B. Fig. ETc) displays Q versus 
/, it shows that a high transmission ratio up to 150 is 
achieved for the working PI channel (i.e. / G [4.0,4.38] 


GW/cm 2 ). The real-time evolution of the optical pulse 
further confirms the high transmission optical diode, see 
Fig. 01 For the case: A B, see Fig. 0(a), not only the 
pulse maintains its soliton-shape after penetrating right- 
CBG region, but also majority of pulse power is trans¬ 
mitted, with the transmissivity up to 53%; while for the 
case: B A, see Fig|4(b), because of the longer CBG 
length, the pulse bounced back, with the transmissivity 
nearly approaching zeros. 

Although pulse dispersion is compensated by the in¬ 
duced optical nonlinearity, resulting in an optical soli¬ 
ton, strong retardation of pulse group velocity is still ob¬ 
served, due to the power conversion between forward- and 
backward-waves in CBG region, see Fig. 2(b). We esti¬ 
mate the soliton propagating velocity from A —>• B, with 
results illustrated in Fig. [3(d). It shows that velocity 
is rather smaller (< 0.03c) for the working PI channel. 
Thus the setting indeed can be termed as a slow-light 
AOSD. 

Furthermore, the breakdown effect of an AOSD based 
on this setting is also investigated. Similar to the break¬ 
down voltage of an electronic diode, when the incident 
intensity I >4.38 GW/cm 2 , the setting becomes two- 
direction conducting. That means such AOSD is broken 
dwon, as indicated from the suddenly drop of Q, see Fig. 
3(c). Fig. [5] demonstrates a typical example of such 
breakdown effect of AOSD. It shows that in both cases 
(A—>B; B—fA), the solitons can penetrate through the 
setting with certain transmission. 

Finally, it should be necessarily pointed out that, in the 
linear case, pulse would experience serious dispersion and 
most of the power is reflected since the pulse wavelength 
is selected near the edge of the bandgap. There is almost 
no propagating wave exhibited inside the gratings. 

Generally, the propagation dynamics of BG solitons 
can be approximately considered as a quasi-particle sys¬ 
tem, which cannot exhibit_ unidirectional transmission 
in an effective potential [31]. However, for the setting 
demonstrated in this work, a Gaussian optical pulse that 
is not an exact solution of the BG soliton cannot develop 
into a BG soliton in CBG segment due to the limited 
CBG length and pulse intensity, see Fig. [2(b). In this 
case, the quasi-particle approximation is invalid in this 
region, and the nonreciprocity of the setting is achieved 
by introducing two spatially asymmetric CBG segments. 
We note that, even though the length of CBG is limited, 
BG solitons can be still generated inside CBG segment 
by launching Gaussian pulses with high-intensity. In this 
case, the quasi-particle approximation is re-built in the 
whole setting, but the diode effect is lost, see Fig. [5j 

In conclusion, we have proposed a scheme to realize 
an AOSD based on a sandwich Bragg-grating structure: 
built as a left-side chirped grating linked to a uniform 
one, and again linked to another different chirped grat¬ 
ing on the right side. Physical parameters were taken for 
silicon. A systematic CMEs has been applied to simu¬ 
late the evolutions of the pulses propagating dynamics 
through the setting. The nonreciprocity of the setting is 
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FIG. 5: (Color online)A typical example of a reverse break¬ 
down for the setting in Fig. [T] The PI of the pulse is I — 4.73 
GW/cm 2 . (a) The pulse injected from side A (n = 0.1 cm), 
propagating through the total setting, (b) The pulse injected 
from side B (7*2 = 0.6 cm), also penetrating the whole setting. 
Here, we fix AS — 7.31 cm -1 and L — 2 cm. 


realized by studying the relationship between the opti¬ 
cal transmission and CBG length. Based on these anal¬ 
ysis, we built an AOSD. A high transmission ratio up 
to 150 was achieved when launching a picosecond Gaus¬ 


sian pulse, with the wavelength located near the edge of 
the photonic bandgap. Because of the balance between 
the optical nonlinearity and pulse dispersion, such pulses 
propagate in a form of solitons, which can preserve its 
soliton-shape during propagating. Moreover, because of 
the introduction of CBG segment, such solitons propa¬ 
gate at a rather small group velocity. Finally, the break¬ 
down effect of AOSD was also investigated. This set¬ 
ting does not require external assistance such as radio¬ 
frequency modulation, magnetic fields, or optical pump¬ 
ing, which make it possible for the photonic applications, 
as well as in the all-optical network and devices produc¬ 
tion, such as optical gate circuit, power isolator, buffer, 
etc. 
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